Trends in Metal Oxide Stability for Nanorods, Nanotubes, and Surfaces 
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The formation energies of nanostructures play an important role in determining their properties, including 
the catalytic activity For the case of 15 different rutile and 8 different perovskite metal oxides, we find that the 
density functional theory (DFT) calculated formation energies of (2,2) nanorods, (3,3) nanotubes, and the (110) 
and (100) surfaces may be described semi-quantitatively by the fraction of metal-oxygen bonds broken and the 
bonding band centers in the bulk metal oxide. 

PACS numbers: 73.22.-f 71.15.Mb 73.21.-b 



The search for cleaner and more sustainable forms of en- 
ergy provides a strong impetus to the development of more 
affordable, active, selective and stable new catalysts to con- 
vert solar radiation into fuels [1,2]. Just as the Haber-Bosch 
process fueled the population explosion of the 20 th century 
[3, 4], it is now hoped that new catalytic processes will pro- 
vide sustainable energy in the 21 st century [1,2]. 

Metal oxides are used extensively as catalysts, electrocata- 
lysts, and photo-electrocatalysts [5-17]. One important prop- 
erty of oxides is their high stability in harsh oxidizing envi- 
ronments compared to their pure metal counterparts. 

The structure of oxide nanoparticles may be determined by 
the surface energy. For instance, for TiC>2 it has recently been 
shown that it is the surface energies which determine whether 
it takes the anatase or rutile structure at the nanoscale [18, 19]. 
Typically, oxide catalysts are in the form of nanoparticles or 
highly porous materials. The catalytic properties of these ma- 
terials are determined to a high degree by the surface [20], and 
control of the surface structure will allow control of the reac- 
tivity [21-24]. However, formation energies of metal oxide 
surfaces are difficult to measure experimentally, and only a 
few values are available in the literature [18, 19]. This makes 
the calculation of trends in surface and nanostructure energies 
an essential first step in understanding the properties of oxide 
catalysts, and the eventual design of novel catalytic materials. 

In this study we provide DFT calculated (2,2) nanorod, 
(3,3) nanotube, and (110) surface formation energies for 15 
different rutile metal oxides, along with the (100) surface for- 
mation energies for 8 different perovskite metal oxides [21]. 
Both the (2,2) nanorod and (3,3) nanotube structures are ob- 
tained by first rolling up a sheet of the material in the PtC>2 
structure. But due to its small diameter, the (2,2) nanorod has 
a bonding structure resembling that of the bulk metal after re- 
laxation, as shown in Fig. 1. 

All DFT calculations have been performed with the plane 
wave code DACAPO using the RPBE exchange-correlation 
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FIG. 1: Structural schematics for (a) Ti02 (2,2) nanorod, (b) TiC>2 
(3,3) nanotube, (c) Ti0 2 (110) surface and (d) SrTi0 3 (100) Ti0 2 
and SrO terminated surfaces. Axes and Miller indices for the rutile 
structures are shown for (a) and (c). 



(xc)-functional [28, 29], converged plane wave cutoffs of 350 
eV and 400 eV for rutile and perovskite metal oxides respec- 
tively, along with a density cutoff of 500 eV. The occupation 
of the Kohn-Sham orbitals was calculated at ksT « 0.1 eV, 
with all energies extrapolated to T = K. The product of 
the supercell dimensions and the number of k-points [30] is 
> 25 A in all repeated directions. For the anti- symmetric per- 
ovskite surface slab calculations a dipole correction has been 
employed [31]. 

We have performed structural relaxations until a maximum 
force below 0.05 eV/A was obtained. At the same time we 
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FIG. 2: Formation energy AEf of rutile metal oxide (2,2) nanorods (•), (3,3) nanotubes (•), and (110) surfaces (■) in eV/MCh vs. (a) bulk 
heat of formation AGj ulk for rutile metal oxides in eV/M02 from Ref. 25, (b) d-band center e& for the bulk metal oxide relative to the Fermi 
energy in eV, and (c) fraction of M-0 bonds broken times the d-band center Ax£d in eV/M02- The experimental nanoparticle surface energy 
[19] and bulk formation energy [26, 27] for rutile T1O2 (□) are provided for comparison. 



have minimized the strain on the unit cell in all periodi- 
cally repeated directions, and employed more than 10 A of 
vacuum between repeated nanorods, nanotubes, and surface 
slabs. Schematics of the T1O2 (2,2) nanorod, (3,3) nanorod, 
and (110) surface, along with the SrTi0 3 (100) Ti0 2 and SrO 
terminated surfaces [32] are shown in Fig. 1. The supercells 
used have been repeated four times along the nanorod and 
nanotube axis in Fig. 1(a) and 1(b), and four times in the sur- 
face plane for the surfaces shown in Fig. 1(c) and 1(d). For 
the TiC>2 (110) surface our four layer thick slab yields surface 
formation energies of 0.44 J/m 2 in reasonable agreement with 
the GGA [33] 12 layer thick slab value of 0.50 J/m 2 [18], and 
the B3LYP [34] 9 layer thick slab value of 0.67 J/m 2 [35]. Dif- 
ferences amongst these values are attributable to the choice of 
xc-functional and number of layers, especially for uneven slab 
calculations, which tend to yield higher energies and converge 
slower [18]. As expected, these values are below the experi- 
mental nanoparticle surface energy of 2.2 ±0.2 J/m 2 for T1O2, 
as shown in Fig. 2 [19]. 

Spin polarized calculations have been performed only for 
CrC>2 and Mn02 rutile metal oxides, while all perovskites 
have been calculated spin polarized. However, spin only 
proved to be important for the energetics of CrC>2, Mn02, 
SrCr0 3 , SrMn0 3 , SrFe0 3 , SrCo0 3 , and SrNi0 3 . A table of 
these results is provided in Ref. 36. 

Serious questions have been raised as to whether standard 
DFT calculatons can describe metal oxides with sufficient ac- 
curacy. For the purposes of the present study, it is worth not- 
ing that it has recently been demonstrated that DFT calcula- 
tions do reproduce semi-quantitatively the formation energies 
of bulk rutile and perovskite metal oxides [25, 37]. 

The formation energies of the various surfaces and nanos- 
tructures from Fig. 1 are shown in Fig. 2. In the following 



we shall analyze the nature of the formation energies, and 
relate them to other characteristics of the materials. An ob- 
vious choice of descriptor for the (110) metal oxide surface 
formation energy Ai^j 110 ^ might be the bulk heat of forma- 
tion AGj ulk , which is well described by DFT for rutile metal 
oxides [25]. However, as shown in Fig. 2(a), these quantities 
appear to be poorly correlated. We instead observe a strong 
correlation between A£^ 110 ^ and the d-band center Ed for the 
bulk rutile metal oxides, as shown in Fig. 2(b). Here, Ed is the 
average energy, relative to the Fermi level ep, of the density 
of states (DOS) projected onto the metal's atomic <i-orbitals in 
the metal oxide rid(e), so that £d = f(s — SF)^d{^)de [38]. 
For the nanorod and nanotube formation energies, we also find 
little correlation with AGj ulk , but a strong correlation with Ed- 

This suggests metal oxide surface and nanostructure for- 
mation energies may be considered perturbations of the bulk 
metal oxide's electronic structure due to bond breaking. In 
this case AEf should be well described by the fraction of 
metal-oxygen (M-O) bonds which are broken A\, times the 
M-0 bond energy in the bulk Em-o- For the case of (110) 
surface formation Ax^ 110 ^ ~ 1 /§ per MO2, while for nanorod 
and nanotube formation A% NT ~ 2 /6 ~ 1 /z per MO2. As- 
suming £m-o ~ —£d, we indeed find that AEf « — A%£d, as 
shown in Fig. 2(c). 

This correlation may be explained qualitatively by recalling 
that the d-band is most stable when the average energy of the 
DOS is £p, i.e. Ed ~ 0. In most cases this is equivalent to 
having a half full d-band [38]. This is nearly the case for 
Ti02, where Ed ~ —0.67 eV, and the (110) surface becomes 
very stable (A£^ 110) w 0.046 eV/Ti0 2 ). As the d-band shifts 
down in energy, the energy stored in the M-0 bonds increases 
accordingly, so that £m-o ~ —£d- 

A more detailed analysis based on the molecular orbitals 
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FIG. 3: DFT calculated formation energy AEf T vs. model prediction AEf odel for (a) (2,2) nanorod (•), (b) (3,3) nanotube (•), and (c) (1 10) 
surface (■) in eV/MCh for rutile and (d) (100) MO2 and SrO terminated surfaces (♦) in eV/SrM03 for perovskites. Open symbols denote 
metal oxides whose most stable phase is not rutile. An error estimate of ±0. 1 eV (grey region) is provided for comparison. 



of the bulk metal oxide is presented in Refs. 36, 39-41. To 
summarize, as the number of d electrons Nd increases, the M- 
O coordinately unsaturated (cus) bond becomes stronger. It is 
only this bond which is broken during (110) surface, nanorod 
and nanotube formation for rutile metal oxides. As shown 
in Ref. 36, this occurs as we progress from non-cus bonding 
(N d =2), M-M bond distortions (N d =3), M-0 cus bonding 
(Nd =4), non-bonding M orbitals with reduced distortions 
(Nd =5,6), to stronger M-0 p a bonds (N d =7). Further, for a 
fixed number of d electrons we also find that the bond strength 
increases as the DOS shifts down in energy. In effect, the 
energy stored in rutile metal oxide M-0 bonds is then equal 
to the energy cost for "shifting" the bonding metal orbitals so 
that their DOS is symmetric about the Fermi level. In such 
a case, both bonding with cus O atoms and distortions of the 
octahedral structure are minimized. 

It is worth noting that we obtain near quantitative agree- 



ment between the DFT calculated formation energies and 
this simple model, when applied to naturally occurring tran- 
sition metal rutile oxides. We now show how this model 
AEy° del = — Ax£d may be extended to post transition metal 
oxides (Ge02, Sn02, and Pb02) and 3d transition metal per- 
ovskites (SrMOa), as shown in Fig. 3. 

For post transition metals, the <i-band is fully occupied, and 
does not participate in the M-0 bonding. Instead, the M- 
O bonding should occur via the occupied portion of the 
metal's p-band. In any case, we find that the M-0 bond en- 
ergy Em-o for post transition metal oxides is approximately 
negative one third the p-band center —e p /3, so that AEf w 
— Axs p /3. This approximation yields near-quantitative for- 
mation energies for post transition metal nanorods, nanotubes, 
and (110) surfaces, as seen in Fig. 3(a), (b), and (c), respec- 
tively. 

Since perovskite metal oxides share the M06 octahedral 
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structure of rutile metal oxides, with interstitial Sr [42, 43], 
it is reasonable to expect that the same correlation with the 
metal d-band centers Sd should hold. From the DFT calcu- 
lations for the perovskite (100) surface, we obtain the aver- 
age of the formation energies for both the MO2 and SrO ter- 
minated surfaces. As such, our model predicts AE^ 00 " 1 « 
-(A X M £™ + A X Sr £ s d r )/2, where A X M ~ 2 / 6 ~ V 3 P er 
SrM0 3 and Ax Sr « 2 / 12 ~ V 6 P er SrM0 3 . From Fig. 3(d), 
we find this is indeed the case. 

In Fig. 3(a), (b), and (c) we also plot the DFT calculated for- 
mation energies versus A£^ odel for transition metals whose 
most stable metal oxide phase is not rutile. For these com- 
pounds, denoted by open symbols in Fig. 3, we find our model 
typically overestimates formation energies. 

For the HexABC layer structure of PtC>2 [21] we find the 
formation energies are completely independent of e& and e p , 
as shown in Ref. 36. However, this is also predicted by our 
model, since no M-0 bonds are broken (A% ~ 0) when form- 
ing the HexABC layer. Instead, the energy cost to form the 
HexABC layer is related to the preference of oxygen to be 
either sp 2 or sp 3 hybridized in the metal oxide. 

In conclusion, we have demonstrated that the formation en- 
ergies for nanotubes, nanorods, and surfaces of metal oxides 
may be determined semi-quantitatively from the fraction of 
M-0 bonds which are broken Ax and the bonding band cen- 
ters £d and e p in the bulk metal oxide. We anticipate such 
models will prove useful in predicting the formation energy 
of doped metal oxide surfaces and nanostructures, their reac- 
tant adsorption energies, and their activities. 
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